Tail States below the Thouless Gap in SNS junctions 
Classical Fluctuations. 



Alessandro Silva 
Department of Physics and Astronomy, Rutgers University, 
136 Frelinghuysen Road, Piscataway 08854, New Jersey, USA. 
(Dated: February 2, 2008) 

We study the tails of the density of states (DOS) in a diffusive superconductor-normal metal- 
superconductor (SNS) junction below the Thouless gap. We show that long-wave fluctuations of 
the concentration of impurities in the normal layer lead to the formation of subgap quasiparticle 
states, and calculate the associated subgap DOS in all effective dimensionalities. We compare the 
resulting tails with those arising from mesoscopic gap fluctuations, and determine the dimensionless 
parameters controlling which contribution dominates the subgap DOS. We observe that the two 
contributions are formally related to each other by a dimensional reduction. 
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I. INTRODUCTION. 

\ The properties of hybrid superconductor-normal metal 
structures (SN) continue to attract considerable atten- 
tion both experimentalljii and theoreticalljiSiSiiiSiSiLSiS, 
though the fundamental process governing the physics of 
such systems, Andreev reflectioniS, has been discovered 
long ago. In fact, while it is well known that generically 
the proximity to a superconductor leads to a modification 

\ of the density of states in the normal metal, the nature 
and extent of this effect depends on the details the hybrid 
structure. In particular, it was recently pointed out^ that 
when a closed mesoscopic metallic region is contacted on 
one side to a superconductor, the resulting DOS turns out 
to depend on its shape. If integrable, the DOS is finite 
everywhere but at the Fermi level, where it vanishes as a 

] power law. On the contrary, in a generic chaotic metal- 
lic region one expects the opening of a gap around the 
Fermi level, the Thouless gap'^. In analogy with the con- 
siderations above, a diffusive metallic region sandwiched 
between two bulk superconducting electrodes has been 
predicted to have a gapped density of states, the gap 
being at energies comparable to to the Thouless energy 
Exh — D/L^, where D is the diffusion constant and Lz 
the width of the normal layeri^i^iSii [see Fig.l]. 

In a diffusive SNS structure with transparent SN in- 
[ terfaces, the density of states in the normal part, aver- 
aged over its thickness, and at energies E right above the 
; gap edge Eg ~ 3A2ETh, is ly (x l/irV ^{E-Eg)/^^ 
' where Aq = {EgS^^^, ^ = l/ii^oV), and V = L^LyL^ 
I is the volume of the normal region. This dependence 
is reminiscent of the density of states at the edge of a 
Wigner semicircle in Random Matrix Theory [RMT] , Aq 
being the effective level spacing right above the gap edge. 
Using this analogy, Vavilov et al.^ realized that the dis- 
order averaged DOS should not display a real gap, but 
have exponentially small tails below the gap edge, anal- 
ogous to the Tracy-Widom tails^^ in RMT. A rigorous 
study in terms of a Supersymmetric Sigma Model de- 
scription of the SNS structure has shown that this is in- 



deed the case^. However, in analogy to the theory of 
Lifshits tails^^ in disordered conductors, the nature of 
the resulting subgap quasiparticle states depends addi- 
tionally on the effective dimensionality d, determined by 
comparing the interface length scales Lx,Ly, with the 
typical length scale of a subgap quasiparticle state, L±. 
In particular, if S> L± > Ly or Lx,Ly ^ L± the 
subgap quasiparticle states are localized either in the x 
direction or in the x — y plane along the interface, respec- 
tively. Correspondingly, the asymptotic tails of the DOS 
deviate from the universal RMT result, applicable only 
in the zero dimensional case [L^, Ly < Lj_]. 

The analogy with RMT applies, within the appropri- 
ate symmetry class, to other physical situations, such 
as diffusive superconductors containing magnetic impuri- 
tie a?i^?d'^ , and superconductors with inhomogeneous cou- 
pling constantsi^. In both cases, at mean field level the 
density of states has a square root singularity close to 
the gap edgeiSiii. Correspondingly, accounting for meso- 
scopic RM-like fluctuation, the disorder averaged density 
of states has tails below the gap edge, with an asymp- 
totics similar to the one calculated in Ref. 9] for SNS 
structures. On the other hand, in the case of diffu- 
sive superconductors containing magnetic impurities, it 
was shown^^'^^ that, in addition to mesoscopic fluctua- 
tions , subgap quasiparticle states can form as a result of 
classical fluctuations , i.e. long-wave fluctuations of the 
concentration of magnetic impurities associated to their 
Poissonian statistics. Similarly, also in superconductors 
with inhomogeneous coupling constant long-wave fluctu- 
ations of the coarse grained gap lead to the appearance 
of subgap quasiparticle states, and consequently to tails 
of the DOS^^. Interestingly, in both cases the tails orig- 
inating from mesoscopic fluctuations and from classical 
ones are formally related by a dimensional reduction^*. 

In this paper, we close this set of analogies, studying 
the contribution to the subgap tails of the DOS in a diffu- 
sive SNS junction arising from long-wave fluctuations of 
the concentration of impurities in the normal layer. Com- 
bining the results of this analysis with those obtained by 
Ostrovsky, Skvortsov, and Feigel'man^i, who considered 
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the subgap tails originating from mesoscopic fluctuations, 
we provide a consistent picture of the physics of the sub- 
gap states. In particular, a quantitative comparison of 
the two contribution shows that mesoscopic fluctuations 
dominate in long and dirty junctions, while classical fluc- 
tuations dominate in wider and/or cleaner ones. In anal- 
ogy with diffusive superconductors with magnetic impu- 
rities, and superconductors with inhomogeneous coupling 
constants, also in the present case the two contributions 
to the subgap tails, arising from mesoscopic and classical 
fluctuations, are related by a dimensional reduction. 

The rest of the paper is organized as follows: in Sec. II 
we present the details of the analysis of the subgap DOS 
arising from fluctuations of the concentration of impuri- 
ties riimp in an SNS junction. In Sec. Ill, we compare the 
two contributions to the subgap DOS associated to meso- 
scopic and classical fluctuations. In Sec. IV, we present 
our conclusions. 



II. SUBGAP DOS ASSOCIATED TO 
FLUCTUATIONS OF ni,„p. 



Let us start considering a diffusive metallic layer in be- 
tween two superconducting bulk electrodes, a geometry 
represented schematically in Fig.l. Assuming kpl » 1, 
where I is the mean free path, this system can be de- 
scribed in terms of the quasiclassical approximation. In 
particular, at mean fleld level [ i.e., neglecting both meso- 
scopic and classical fluctuations ], neglecting electron- 
electron interaction, and assuming the thickness of the 
metallic layer >> I, one can describe the SNS struc- 
ture by the Usadel equationSSiSi 



^V^0 + i E sm\9] = 0, 



(1) 



where D = VpT/3 is the diffusion constant, E is the 
energy measured from the Fermi level, assumed to be 
I £^ 1^ A, where A is the gap in the bulk electrodes. 
The field 9 is related to the quasiclassical Green's func- 
tions and the anomalous Green's function by the relations 
g{r,E) = cos[6'(r,£;)], fir,E) = ism[e{r, E)]. In addi- 
tion, assuming the interfaces to be perfectly transparent, 
the proximity to the two superconducting regions can be 
described by the boundary conditions 9{z — ±L^/2) = 
tt/2. 

It is convenient to measure all lengths in units , and 
set 9 = tt/2 + i'i' . Therefore, Eq.QJ becomes 



E 
) 

Exh 



cosh[*] = 0, 



(2) 



where Eth = D/L1 is the Thouless energy. The bound- 
ary conditions for the fleld are simply '^{z — ±1/2) = 
0. 

In terms of 4* the DOS is = 2i/oIm[sinh[\E']], where vq 
is the density of states of the normal metal at the Fermi 
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FIG. 1: A schematic plot of an SNS junction: two bulk super- 
conducting electrodes (S) connected to a diffusive metal (N) 
of thickness L,. The interfaces have linear size Lx^ Ly. 



level. The DOS can be calculated looking for solutions 
of Eq.||2Jl uniform in the x — y plan o"^i^i^i^ . In particular, 
for E <Eg = C2ETh [C2 ^ 3.122] aU solutions of Eq.Q 
are real, implying = 0. Therefore, one identifies Eg 
with the proximity induced gap within the normal metal 
layer. The mean field DOS right above Eg averaged over 
the z direction is found to be 



3.72 1^01 



IE-E„ 



E„ 



(3) 



Let us proceed analyzing the tails of the DOS at en- 
ergies E < Eg arising from fluctuations of the con- 
centration of impurities, i.e. long-wave inhomogeneities 
in the x — y plane of 1/r. We first consider an SNS 
structure such that the linear size of the SN interfaces 
is much larger than the thickness of the metallic layer 
[Lx,Ly ^ Lz]- In the framework of the Usadel de- 
scription of the metallic layer [Eq.jSJ] one can account 
for long-wave transversal fiuctuations of the concentra- 
tion of impurities by promoting Exh, or equivalently 
Eg = C2ETh, to be a position dependent random vari- 
able, characterized by the statistics 



(SEgi^)) 
{5Eg{^)6Eg{^')) 



Eg+5Eg{jC), 

0, 



^J(x-x') 



(4) 
(5) 

(6) 



where d is the effective dimensionality of the system, and 
Ud the effective concentration of impurities. As shown be- 
low, d is determined by comparing the linear sizes of the 
interface Lx,Ly to the linear scale of the subgap states 
~ Lj{{Eg - E)lEgfl^. If L.^,Ly > the system 
is effectively two dimensional, and ~ nimpLz- On the 
other hand, if < L± ^ Ly [or Ly < L± <C L^], 
the system is effectively one dimensional, and ni — 

Accounting for these fiuctuations, the Usadel equation 
Eq.(|2l becomes 

dl^ + V^vl/ + 2C2^(1 - <5e<,(x)) cosh[*] = 0, (7) 
where Stg = 5Eg/Eg. 
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Our purpose is to calculate the DOS averaged over fluc- 
tuations of 6Eg at energies E < Eg. For this sake, let us 
introduce SE = Eg- E, and S^iz, x) = x) - *o(^), 
where is the solution of Eq.® a.t E — Eg. Expanding 
Eq.lUI) and keeping the lowest order nonlinearity in i5\l/ 
one obtains 



— — yi)y^jyu<^- u^gj 



go{z){Se-Seg),i8) 



where Se — SE/Eg, ga{z) = 2C2 cosh[^'o(2:)], and fa{z) = 

In order to simplify further Eq. JH)) , it is useful to notice 
that the operator H = —d^ — /o(^), diagonalized with 
zero boundary conditions at ±1/2, admits an eigenstate 
$0 with zero eigenvalue. Physically, $0 determines the 
shape of the mean field z-dependent DOS obtained from 
Eq.(|2J). Therefore, it is natural to set 



J*(z,x) ~ y/Ai/A2 x(x) $o(^), 



(9) 



with Ai = J dz go $0 ~ 7.18, and A2 = / ^ $g ~ 
2.74. 

Substituting Eq.© in Eq.JSJ, and projecting the re- 
sulting equation on $01 one obtains 



V\ + ^ Se - Segi^) 



(10) 



where we rescaled the length by {Ai A2) and 

(56,(x)<5e,(x')) =?7<5(x-x'), (11) 

withr,= (AiA2)i/V(nd Li). 

Let us now split x = ~u + ^-iid obtain the system 



-V^u + - = 6e - 6e„, 



-V^ V + UV ^0. 



(12) 
(13) 



Interestingly, this set of equations is analogous to the 
equations obtained by Larkin and Ovchinikov in the con- 
text of the study of gap smearing in inhomogeneous su- 
perconductorsii, and to the equations obtained by the 
author and loffe in the context of the study of subgap 
tails in diffusive superconductors containing magnetic 
impurities 

Let us now proceed with the calculation of the DOS. In 
the present notation, the DOS averaged over the thick- 
ness of the normal layer is given by 



i^(x, (5e I Seg{-x.)) 



3.72w(x,<5e I 5eg(x)). (14) 



We are interested in calculating the average density of 
states — 3.72{v) at energies below the Thouless 

gap [Se > 0]. In this parameter range, the corresponding 
functional integral 

^ / D[SegH^,5e \ j6,(x))e-^/(^") '^''^-''^"^ ^(^^^ 



receives its most important contributions by exponen- 
tially rare instanton configurations of deg such that, at 
specific locations along the interfaces of the junction, 
'^£g(x) > Se. The remaining task is to select among all 
these fiuctuations the one that dominates the functional 
integral Ea. (|15|l . i.e. the optimal fluctuation . 
The action associated to a configuration of 5eg is 



u-u^ + 6ef, (16) 



where we used Ea. (|12|l to express Seg in terms of 
and, with exponential accuracy, neglected the term 
in the action. In order to find the optimal fluctuation 
one has to find a nontrivial saddle point uq of S, tending 
asymptotically to the solution of the homogeneous prob- 
lem [uq — > VSe], and subject to the constraint of having 
nontrivial solutions for v of Eq. (|13|l . 

Since the normal metal layer is diffusive, and momen- 
tum scattering isotropic, it is natural to assume the opti- 
mal fiuctuation to be spherically symmetric. The Euler- 
Lagrange equation associated to S is 



1 



where 



(_ A^'^) -Ku) (A^'') u — u + Se) — 



(17) 



(18) 



is the radial part of the Laplacian in spherical coordi- 
nates. An obvious solution to Eq. H17() is obtained setting 



A^'^^u-u^ + 5e = 0. 



(19) 



This equation is equivalent to the homogeneous Usadel 
equation with uniform Eg. i.e. Ea. (|10|l with Seg = 0. 
Though this equation has definitely nontrivial instanton 
solutions for u with the appropriate asymptotics, it is 
possible to show that the constraint of Eq. (|13l) is satisfied 
only by f — 0. This is physically obvious since Ea. l|19() 
describes a uniform system where all long-wave fluctu- 
ations of 1/t have been suppressed, and thus, within 
the present approximation scheme, the subgap DOS must 
vanish. However, it should be pointed out that, account- 
ing for mesoscopic fluctuations, the instanton solutions 
of Eq. lfTO|) describe the optimal fluctuation associated to 
mesoscopic gap fluctuations, as shown in Ref.0. 

Let us now look for the nontrivial saddle point. Equa- 
tion (|17|) is equivalent to the system 



(-lA^'i)+u)h = 0. 
A^'^'>u-u^ +Se^ h. 



(20) 
(21) 



which can be reduced to a single second order instanton 
equation setting h = {2dru)/r. With this substitution, 
Ea. (|20|) becomes the derivative of Ea. H21|) . which now 
reads 



Se^O. 



(22) 
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Notice that this equation is, upon reduction of the dimen- 
sionaUty by 2, identical in form to the one associated to 
mesoscopic fluctuations, Ea. lfH^ . As we will see later, 
this reduction of dimensionality relates in a similar way 
the dependence of the action associated to classical and 
mesoscopic fluctuations on 5e. 

It is now straightforward to see that the instanton so- 
lution uq of this equation with the appropriate asymp- 
totics describes indeed the optimal fluctuation, the con- 
straint of Eq. (|13|l being automatically satisfied in virtue 
of Ea. H2U|l . with vq {2drUQ)/r. Moreover, the corre- 
sponding optimal fluctuation of Seg is deg = 2drUf)/r. 

It is clear that the instanton solutions of Eq. (|22ll must 
have the form uq = \f5ev{r / X), with A = l/((5e)^/'*. The 
corresponding equation for v{r) is dfv + {d — 3)/rdrV — 
v'^ + l — 0. The instanton solution of this equation can be 
easily found numerically, and the corresponding action S 
calculated. The result is 



Sd = CLdiidLt 5e 



(23) 



where the constants ad are oi ~ 0.88, and 02 ~ 7.74. 

Within our approximation scheme, the density of 
states is {v) oc W ex-p[—S], where is a prefactor due to 
gaussian fluctuations around the instanton saddle point. 
The calculation of W can be performed using the stan- 
dard technique due to Zittarz and Langer, and is similar 
to those reported in Ref. 18 17J. To leading order in the 
saddle point approximation, the final result is 



1^0 



Pd \ rid if 5t 



d(lQ--d)-12 



(24) 



where /3i « 0.1 and j32 ~ 0.5. 

The result in Eg. 1)24(1 relies on a saddle point approx- 
imation, which is justified provided 5^ S> 1. This trans- 
lates into the condition 



5e > 



1 



MTidL'i 



(25) 



As mentioned before, the effective dimensionality, and 
therefore the asymptotic density of states, is determined 
by comparing the linear size of the optimal fluctuation, in 
dimensionfull units — L^X — Lz/Se^/^, to the linear 
dimensions of the interfaces L^^ Ly. If L^, Ly ^ L± the 
asymptotics is effectively two dimensional [d = 2], while 
for Ly ^ L±, L± the asymptotic DOS is effectively 

one dimensional [d = 1]. Since L± increases as the energy 
gets closer to the average gap edge, it is clear that in 
any flnite size system the applicable asymptotics might 
exhibit various crossovers, 2d — > Id ^ Od, as Se 0. In 
particular, the tails are zero dimensional when L^, Ly < 
L±, in which case the asymptotic form of the DOS is 
obtained by calculating the integral 



1^0 



^ ~ 3.72 



djSeg) 
\/27r?7o 



\/Seg — Se e 



-So 



(26) 



where rjo 



LxLyLz] and 5*0 



III. MESOSCOPIC VS. CLASSICAL 
FLUCTUATIONS. 

In the previous section we have discussed the asymp- 
totic density of states below the Thouless gap originating 
from classical fluctuations, i.e. inhomogeneities in the 
concentration of impurities or equivalently in l/r. As 
discussed in the introduction, this mechanism to gener- 
ate subgap states is complementary to mesoscopic fluc- 
tuations of the gap edge. 

The tails associated to mesoscopic gap fluctuations 
have been calculated by Ostrovsky, Feigel'man and 
Skvortsov in Rcf. 9] . To exponential accuracy, the subgap 
DOS associated to mesoscopic fluctuations is (i^) /vo oc 
exp[— S'd], where 



Sd ^ dd Gd (Se)' 



(27) 



where dd is a constant [do ~ 1.9, di ~ 4.7, and 02 — 10], 
and Gd is the effective dimensionless conductance 



Go = iTTj^oD^.^ 

Gi = AttvqDLx, 
G2 = AttvqDLz. 



(28) 

(29) 
(30) 

The scale of the optimal fluctuation associated to meso- 
scopic fluctuations is also L_l ~ L^/ {SeY^'^ ■ Therefore, 
the effective dimensionality d is to be determined accord- 
ing to the criteria presented in the previous section. 

Before discussing the comparison of mesoscopic and 
classical fluctuations, let us flrst explain the rationale 
behind the separation these two contributions. Though 
it is clear that the only physical fluctuations in a real 
sample are associated to fluctuations in the positions of 
impurities, these fluctuations can affect the DOS in two 
ways: i)- depress the Thouless gap edge by increasing 
locally the scattering rate [classical fluctuations], or ii)- 
take advantage of interference effects in the quasiparti- 
cle wave functions to generate quasiparticle states that 
couple inefficiently to the superconducting banks [meso- 
scopic fluctuations] . It makes sense to think of two types 
of effects separately if the actions associated to them are 
very different in magnitude [S* ^ 5 or vice versa]. Ob- 
viously, in the crossover region, where S ~ S the sepa- 
ration of these two mechanisms is meaningless, because 
the system can take advantage of both at the same time. 

With this caveat, let us proceed in the comparison 
of these two contributions, starting with the zero di- 
mensional case. Since the dimensionless conductance is 
Go « Eg/6, where 6 « l/{voV) is the level spacing, then 
the d = action associated to mesoscopic fluctuations 
can be written as 



3/2 



(31) 
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where Aq = {EgS"^)^^^, where 5 = 1/(1/0 V) is the level 
spacing in the metallic layer. Physically, Ap can be in- 
terpreted as being the effective level spacing right above 
the gap edge. Indeed, from Eq.© one sees that 



(32) 



Therefore, the result of Ea. (|31|) indicates that tails orig- 
inating from mesoscopic fluctuations of the gap edge are 
universal [in d = 0], in accordance to the conjecture for- 
mulated in Ref . ^ on the basis of Random Matrix Theory. 
In turn, in the zero dimensional case the action associ- 
ated to classical fluctuations is 




^0 



SE 
SE^ 



(33) 



where SEq = Eg/ y/riimpV is the scale of typical fluc- 
tuations of the gap edge associated to fluctuations of 
the concentration of impurities. The dimensionless pa- 
rameter controlling which which mechanism dominates 
is therefore 



70 = 



A^ 

SEo' 



(34) 



Clearly, for 70 ^ 1 mesoscopic fluctuations dominate the 
subgap tails, while for 70 ^ 1 classical fluctuations give 
the largest contribution to the subgap DOS^"*. 

Let us now write 70 in terms of elementary length 
scales, one can estimate 



1 



70 



7/6 



1 

M (L^Ly/py/e 

kpl (L^Ly/Pf/^' 



(35) 



where we used the fact that the scattering cross section 
of a single impurity a is typically of the same order of 
A|^. Within the assumptions of the theory, 70 is the ratio 
of two large numbers, and therefore its precise value de- 
pends on the system parameters. However, from Ea. H35|) 
we see that making the junction longer and longer, i.e. 
increasing ix, tends to favor mesoscopic fluctuations. In- 
tuitively, this is due to the fact that as increases, 
the dimensionless conductance of the junction diminishes 
while the average number of impurities increases, there- 
fore suppressing the associated fluctuations of the gap 
edge. At the same time, increasing the area of the junc- 
tion, or making them cleaner, reverses the situation. In 
summary, mesoscopic fluctuations are favored in long and 
dirty junctions, while classical fluctuations are favored in 
wider and/or cleaner ones. 

Since in higher dimensionalities the linear scale of the 
optimal fluctuation associated to the two mechanism is 
identical [L±^ = Lz/{Sey^^], it is possible, and physically 
suggestive, to reduce the form of the actions in d = 1, 2 



to a zero dimensional action calculated within the typical 
volume of the optimal fluctuation. The latter is V± — 
LxLj^Lz for d = 1, and V±_ = L^L^ m d = 2. For 
example, for d — 1 one can write 



where SE. 



Si 



^E, 



SE 



SE, 



''imp 



V±. Similarly, 



^1 



SE 



(36) 



(37) 



where Ae// = {EgSlffY/^, Seff = 1/{i^qV±) being the 
level spacing in the volume of the optimal fluctuation. 
In analogy to the zero dimensional case, one is therefore 
led to conclude that also in for one dimensional tails long 
and dirty junctions are dominated by mesoscopic fluctu- 
ations, while wider and/or cZeaner junctions favor clas- 
sical ones. This qualitative statement is indeed correct, 
but the proof is complicated by the energy dependence 
L^. 

The appropriate way to proceed for d ~ 1, 2 is to write 
the actions associated to classical and mesoscopic fluctu- 
ations in compact form as 



(38) 
(39) 




i?<,/(arfn<iLf)4/(8-'i), and A^ 



where SEd — 
Eg/{ddGd)'^^^^^'^^ ■ Therefore, the dimensionless param- 
eter that determines which contributions dominates the 
subgap DOS is 



Id 



Ad_ 

SEd 



(40) 



If 7d ^ 1, the subgap DOS is dominated by mesoscopic 
gap fluctuations, and the applicable result is Ea. (|27|l . On 
the other hand, for 7d ^ 1 the DOS below the gap is de- 
termined by long-wave fluctuations of 1/t [Eq.(|2i|l]. Fi- 
nally, estimating 7^ in terms of elementary length scales, 
one obtains 



71 
72 



1 {Lz/l)'/' 



(41) 

(42) 



In analogy to Eq. H35() , the fact that ^d is proportional to 
a power of Lz/l implies that mesoscopic fluctuations are 
dominant in long junctions, while the inverse proportion- 
ality of 7d on a power oik pi and of the linear size of the 
interface [in c? = 0, 1] implies that wide interfaces and/or 
cleaner samples may favor the contribution arising from 
classical fluctuations. 
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IV. CONCLUSIONS 

In this paper, we discussed the effect of inhomogeneous 
fluctuations of the concentration of impurities, or equiv- 
alently of 1/r, on the tails of the DOS below the Thou- 
less gap in diffusive SNS junctions. We have shown that 
these classical fluctuations lead to the formation of sub- 
gap quasiparticle states and are complementary to meso- 
scopic fluctuations in determining the asymptotic DOS. 
Finding the dimensionless parameter that controls which 
mechanism gives the dominant contribution to the sub- 
gap tails, one finds that, qualitatively, mesoscopic fluc- 
tuations are favored in long and dirty junctions, while 
classical ones dominate in wider and/or cleaner ones. 

We have observed that, as for diffusive superconduc- 
tors containing magnetic impurities, and for diffusive su- 
perconductors with an inhomogeneous coupling constant. 



the two contributions are formally related by a dimen- 
sional reduction by 2, both at the level of instanton equa- 
tions determining the optimal fluctuation, and in the de- 
pendence of the DOS on the distance from the gap edge 
6e. As in other physical systems^'^, it is natural to expect 
that supersymmetry is at the root of dimensional reduc- 
tion also in this context. This fact could in principle 
be elucidated generalizing the Sigma Model describing 
mesoscopic fluctuations to include the physics associated 
to classical fluctuations. 
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